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LONG-TERM  GOALS 

The  long-term  objective  of  our  research  is  to  quantify  the  structure  of  turbulence  in  fluvial  and 
estuarine  environments,  in  order  to  develop  remote-sensing  tools  for  environmental  assessment  as  well 
as  to  improve  numerical  simulations. 

OBJECTIVES 

The  objectives  of  this  program  were: 

•  to  quantify  the  turbulence  length  scale  and  turbulent  dynamics  in  an  estuary  under  varying 
stratification  conditions  and  geometries,  including  relatively  uniform  boundary-layer  flows  and 
highly  disrupted  wake  flow  conditions; 

•  to  quantify  the  key  properties  of  observed  coherent  structures,  including  horizontal  and  vertical 
scales,  intensity  of  vertical  motions, 

•  to  provide  a  field-scale  test  of  turbulence  closure  models  and  large-eddy  simulations  via  direct 
measures  of  turbulent  kinetic  energy,  length  scale  and  turbulent  dissipation  rate  combined  with 
accurate  measures  of  the  Reynolds-averaged  quantities. 

•  to  work  with  the  other  COHSTREX  investigators  to  ascertain  the  relationship  between  surface 
signatures  of  the  turbulence  and  the  characteristics  of  the  Reynolds-averaged  flow,  turbulence, 
density  structure  and  bathymetry. 

APPROACH 

The  Mobile  Array  for  Sensing  Turbulence  (MAST)  is  a  10-m  long  aluminum  structure  with  6  sets  of 
instruments  for  measuring  the  turbulent  quantities  and  Reynolds-averaged  velocity  and  density  at 
multiple  depths.  The  instrument  is  deployed  off  the  side  of  a  research  vessel  (fig.  1 ),  either  in 
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Figure  1.  The  MAST  turbulence  measurement  system,  a)  in  stowed  configuration  on  the 
Centennial  in  the  Snohomish  River  (June  2006);  b)  schematic  of  deployed  configuration;  c) 
up  of  top  sensor  bracket;  d)  schematic  of  sensors  in  deployed  configuration. 
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underway  mode  or  at  anchor.  The  depths  of  the  sensors  are  varied  by  the  tilt  of  the  mast,  with  the 
sensor  spacing  1-1.5  m  in  the  vertical.  Turbulence  quantities  are  measured  with  co-located  Seabird 
SBE-7  micro-conductivity  sensors  and  Sontek  ADVs  (fig.  1).  The  stratification  is  measured  with 
RBR  T-S  sensors.  The  spatial  resolution  of  turbulence  measurements  is  5-10  cm  for  velocity  (limited 
by  the  25  Hz  sampling  rate  of  the  ADVs  and  flows  past  the  sensors  of  50-100  cm/s),  and  less  than  1  cm 
for  conductivity  (effective  sampling  rate  200  Hz). 

The  first  field  deployment  of  the  MAST  was  the  2006  COHSTREX  experiment  in  the  Snohomish 
River  (Fig.  2).  There  were  three  modes  of  observation:  one  was  a  set  of  anchored  measurements  in 
close  proximity  to  a  submerged  jetty  where  Jessup  and  colleagues  collected  infrared  observations;  the 
second  was  a  30-hour  anchor-station  to  observe  the  transitions  between  unstratified  and  stratified 
boundary-layer  turbulence  in  the  estuarine  channel,  and  the  third  mode  was  a  set  of  along-river  and 
across-river  transects  of  the  velocity  and  water  properties. 
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Figure  2.  The  Snohomish  River  estuary,  showing  the  locations  of  the  anchor-station 
observations.  The  triangle  indicates  a  uniform  reach  of  the  estuary,  where  stratified  and 
unstratified  boundary-layer  conditions  were  observed.  The  circles  indicate  observations  near  the 
jetty,  where  “wake”  conditions  were  observed. 

The  turbulence  analysis  follows  Kaimal  et  al.  (1972),  fitting  the  velocity  and  conductivity  spectra  to 
universal  shape  functions,  in  order  to  identify  the  turbulent  length  scale  and  to  quantify  the  turbulent 
kinetic  energy  (tke)  and  dissipation  rate.  The  vertical  velocity  spectrum  was  found  to  provide  the 
cleanest  fits  to  the  Kaimal  spectrum,  and  most  of  the  length  scale  analysis  focused  on  the  vertical 
velocity  spectrum.  With  reliable  estimates  of  the  turbulent  length  scale,  the  turbulent  kinetic  energy 
and  the  dissipation  rate,  the  scaling  of  turbulence  for  different  flow  conditions — unstratified  boundary 
layers,  stratified  boundary  layers,  and  wakes — could  be  tested. 

RESULTS 


Methodological  Advances 

The  study  confirmed  that  the  MAST  provides  an  effective  tool  for  quantifying  turbulence  in  estuarine 
and  fluvial  environments.  The  inertial  subrange  method  is  well  suited  to  the  high-energy  conditions  of 
rivers  and  estuaries,  due  to  the  high  Reynolds  number  of  the  turbulence  and  the  correspondingly  wide 
intertial  subrange.  The  MAST  is  well  suited  to  quantifying  vertical  distributions  of  the  Reynolds 
averaged  and  turbulence  quantities,  which  is  essential  for  determining  the  influence  of  stratification 
and  proximity  to  boundaries.  It  is  also  effective  for  quantifying  buoyancy  and  momentum  fluxes, 
which  in  conjunction  with  the  measurements  of  Reynolds-averaged  vertical  gradients  provide  explicit 
tests  of  turbulence  closure  formulations. 

Turbulence  Length-Scale 

The  Kaimal  spectral  analysis  provides  robust  estimates  of  the  turbulence  length-scale  Lw,  based  on  the 
horizontal  wavenumber  spectrum  of  vertical  velocity  variations  (hence  the  subscript  w;  Fig.  3).  The 


Example  Velocity  Spectrum 


Figure  2.  Wavenumber  spectra  of  vertical  (blue)  and  horizontal  (green)  velocity  during  the 
Snohomish  field  study.  The  Kaimal  fit  to  the  vertical  velocity  spectrum  is  shown,  with  the  roll¬ 
off  wavenumber  ko  indicated.  The  turbulence  lengthscale  is  estimated  as  Lt=  0.56/k„. 


length-scale  should  provide  a  key  diagnostic  as  to  the  nature  of  the  turbulence.  Specifically,  U 
should  vary  as  the  distance  from  the  boundary  for  unstratified,  boundary-layer  turbulence,  and  it 
should  vary  as  the  Ozmidov  scale  Lo=  (e/N^)''^^  within  stratified  shear  layers.  Deviations  from  these 
scaling  relationships  should  indicate  non-equilibrium  conditions,  for  example  wakes  and  fronts. 

The  length-scale  analysis  confirmed  that  in  “normal”  circumstances,  i.e.,  relatively  homogeneous  and 
stationary  flow  conditions,  the  turbulence  length-scale  is  set  by  either  the  Ozmidov  scale  or  the 
distance  to  the  boundary,  depending  on  which  is  smaller.  In  figure  3,  the  Ozmidov  scale  is  plotted 
against  the  MAST  estimate  of  the  turbulent  length  scale,  where  both  quantities  have  been  non- 
dimensionalized  by  boundary  layer  scaling.  These  two  quantities  show  roughly  a  linear  relationship  for 
values  on  the  x-axis  less  than  one,  where  the  Ozmidov  scale  is  smaller  than  boundary  layer  scaling. 

The  data  asymptote  to  a  value  of  roughly  one  where  the  Ozmidov  scale  exceeds  the  distance  to  the 
boundary,  consistent  with  traditional  boundary  layer  theory.  This  result  provides  a  confirmation  that 
the  MAST  can  reliably  quantify  the  turbulent  length  scale,  and  it  also  provides  field  confirmation  for 
the  Ozmidov  scaling  of  the  length  scale  in  highly  stratified  conditions,  for  which  there  are  few  field 
observations. 
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Figure  3.  Turbulence  lengh-scale  non-dimensioanalized  against  the  boundary-layer  length-scale 
Lbi  =Kz(l-z/h)  (y  axis),  compared  to  Lo  /Lbi  (x-axis),  for  the  uniform  boundary  layer  (a)  and  the 
wake  of  the  jetty  (b).  The  uniform  boundary  layer  data  confirm  the  conventional  scaling  for  the 
lenthscale,  whereas  the  wake  data  indicate  that  the  turbulence  lengthscale  is  generally  greater  by 

roughly  a  factor  of  2. 

The  data  collected  in  the  vicinity  of  the  jetty  provide  a  marked  contrast  to  the  data  collected  during  30- 
hour  anchor  station.  MAST  estimates  of  the  turbulent  length  scales  exceed  those  expected  based  on 
boundary  layer  theory  by  roughly  a  factor  of  2  (figure  3).  Near  surface  estimates  of  length  scale  are  in 
good  agreement  with  estimates  obtained  from  the  infrared  observations  by  Jessup  and  colleagues. 
Whereas  the  30-hour  anchor  station  data  are  consistent  with  either  Ozmidov  or  wall  layer  scaling,  the 
presence  of  the  jetty  imposes  an  additional  length  scale  on  the  flow,  increasing  the  turbulent  length 
scale  relative  to  the  local  production  scale. 

Another  departure  from  conventional  scaling  for  the  lengthscale  was  observed  during  the  passage  of 
the  salinity  front  during  the  flooding  tide.  The  turbulence  scale  exceeded  the  Ozmidov  scale  by 
roughly  a  factor  of  two,  indicating  non-equilibrium  turbulence  conditions. 

Both  the  jetty  data  and  the  observations  in  the  front  indicate  that  in  non-equilibrium  conditions,  the 
turbulence  lengthscale  may  significantly  exceed  its  equilibrium  value.  Both  the  jetty  data  and  the 
observations  in  the  front  indicate  that  in  non-equilibrium  conditions,  the  turbulence  lengthscale  may 
significantly  exceed  its  equilibrium  value.  These  results  are  consistent  with  the  hypothesis  that 
turbulence  production  exceeds  dissipation  in  regions  of  rapid  transitions  in  flow  conditions,  either  as  a 
result  of  changes  in  flow  geometry  or  the  presence  of  fronts. 

IMPACT/APPLICATIONS 

The  MAST  provides  unprecedented  resolution  of  turbulence  length  scale  and  turbulent  kinetic  energy 
through  the  water  column  in  a  coastal  environment.  Whereas  microstructure  measurements  (e.g., 

Peters  and  Bokhorst,  2000)  are  effective  for  estimating  dissipation,  they  provide  little  information 
about  the  energy-containing  scales  of  the  turbulence  or  the  total  tke.  The  ability  of  the  MAST  to 
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obtain  continuous  timeseries  at  fixed  levels  provides  the  statistical  information  required  to  obtain  the 
information  about  the  outer  scales  of  the  turbulence. 

The  key  seientific  results  of  this  study  are  1)  eonfirmation  of  the  Ozmidov  scaling  of  the  turbulence 
lengthscale  for  uniform,  stratified-flow  conditions,  but  2)  significant  deviations  from  either  boundary- 
layer  or  Ozmidov  scaling  in  regions  of  flow  transitions,  eonsistent  with  local  excess  of  turbulenee 
production.  These  observations  indicate  that  non-equilibrium  turbulence  is  important  at  the  scales 
relevant  to  rivers  and  estuaries,  and  therefore  turbulence  closure  must  properly  account  for  departures 
from  a  production-dissipation  balance.  Second-order  turbulence  closure  models  include  non¬ 
equilibrium  dynamics,  but  few  measurements  have  heretofore  been  available  to  test  model 
performance  under  non-equilibrium  conditions.  Additional  analysis  of  these  data  will  provide  a 
valuable  test  of  second-order  closure  models. 
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